A Scalar Wigner Theory for Polarized Light in Nonlinear Kerr Media 
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A scalar Wigner distribution function for describing polarized light is proposed in analogy with the treatment 
of spin variables in quantum kinetic theory. The formalism is applied to the propagation of circularly polarized 
light in nonlinear Kerr media and an extended phase space evolution equation is derived along with invariant 
quantities. We further consider modulation instability as well as the extension to partially coherent fields. 
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I. INTRODUCTION 

The field of statistical optics, and its corresponding quan- 
tum statistical version, is a field where basic science ques- 
tions meet applications, the latter in the form of, e.g., optical 
communication systems. Formalisms developed for the quan- 
tum realm has been transferred to classical problems, and the 
strong coupling to kinetic descriptions of plasmas and quan- 
tum kinetics of particles has given rise to interesting cross- 
breeding between seemingly disparate areas of research [ 1 - 
3]. The description of atomic, molecular and optical physics, 
in particular quantum optical systems, in terms of statistical 
quantities lends itself to a multitude of applications, due to 
its direct relevance for the interpretation of experimental re- 
sults (see, e.g., [4]). The goal here is less ambitious; we are 
interested in the possibility of formulating and interpreting a 
scalar (quasi-)distribution function of polarized light. The ex- 
tension of our results to applications will be left for future 
studies (something we discuss further in the Conclusions). 

The spin and orbital angular momentum of light has over 
the last few years attracted an increasing interest, for a num- 
ber of reasons. One such aspect lies in the possibility to extract 
more information from light sources, taking the angular mo- 
mentum of light into account J5|] . There has also been interest 
in this field due to the possibility to manipulate material sys- 
tems through optical tweezers centered around tailored light 
sources with angular momentum @]. Thus, the dynamics of 
the polarization state of light, and its interaction with mat- 
ter, lies at the heart of some recent applications. Moreover, 
the spin degrees-of-freedom is something that recently has at- 
tracted interest in the field of quantum plasmas If7l-fl0ll. much 
due to the fact that it shows a potential to lead to new and im- 
portant applications (see, e.g., (HHMl). We will here adopt 
an approach common to some of these works, in particular the 
approach of obtaining a scalar distribution function from vec- 
tor degrees-of-freedom using a suitable projection operator. 
Since we below will deal with what is essentially a two-state 
system, the natural projection operator follows closely from 
the spin-1 /2 case. 

The setup in this paper is based on the nonlinear interaction 
of light, through a Kerr medium. Thus, the founding descrip- 
tion is given by a coupled set of nonlinear Schrodinger equa- 
tions (NLSE). The NLSE has an enormous breadth of appli- 
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cation, and can be found essentially whenever a wave system 
experiences slow modulations lfl5lll6ll . Due to its structure, it 
naturally exhibits soliton formation, self-focusing, and wave 
collapse lfl5ll . Thus, it can be used for describing optical lfl7ll 
and ocean rogue waves IU8I1 . the propagation of light through 
a nonlinear quantum vacuum lfl9l l20ll. solitary waves in opti- 
cal fibers ll2lll . atmospheric waves l22ll . and the formation of 
structures in plasma waves 12311 . Thus, the NLSE allows for 
a very encompassing formulation which can be applied to a 
large variety of fields. 

In this paper, we investigate the statistical dynamics of po- 
larized light in a nonlinear refractive medium, starting from 
a two-state system described in terms of two coupled non- 
linear Schrodinger equations (NLSEs). In particular, we ob- 
tain a generalized Wigner-Moyal type evolution equation for a 
scalar quasi-distribution function, where the polarization state 
is described by an extra independent variable. Here differ- 
ent polarization states corresponds to distinct positions on the 
Poincare sphere for this variable. The conservation laws for 
our scalar distribution function are analyzed and the Vlasov 
limit is described. Furthermore, the modulational instability is 
studied for a fully coherent case, and a simple generalization 
to partially coherence is considered. Finally the applications 
and possible generalizations of our results are discussed. 

II. FUNDAMENTALS 

Two circularly polarized modes can be used to form a ba- 
sis for a general polarization state. In this paper we consider 
propagation of two counter-rotating modes of right- and left- 
handed circularly polarized light in a guiding nonlinear Kerr 
medium, such as e.g. an optical fiber. The two modes are 
assumed to obey the following coupled system of nonlinear 
Schrodinger type evolution equations, c.f. lf24T,l2~5ll . 
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where we have assumed the propagation direction to be along 
the z-axis with x representing either the transverse coordinate 
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or time depending on the situation under study. Correspond- 
ingly a represents either the diffraction or dispersion coef- 
ficient while v g ± represents the group-velocity of the corre- 
sponding mode and y determines the strength of the nonlin- 
earity. The system further includes both linear and nonlinear 
coupling with coupling coefficients K and v, respectively. For 
a low birefringence medium we have j3+ = j3_ = j5 x + j3 v and 



K = f3 x — j8y, with f} x (fly) being the linear propagation constant 
of the slow (fast) axis. 

The system Eqs.ffl} and (O has been extensively studied 
in the context of nonlinear pulse propagation in optical fibers 
ll25ll . It is a conservative system and can be shown to have the 
following Lagrangian 
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The aim of this paper is to introduce a scalar Wigner func- 
tion and a concomitant evolution equation equivalent to the 
system Eqs.(Q]) and {2). To this end we define a correlation 
matrix T mn (x 1 , x 2 , z) = ( y* m {x\ , z) \jf n (x 2 , z) ) and a Wigner ma- 
trix W mn = W mn (x,p,z), viz. 

W,„„(x,p,z) = ^J T mn {x + ^/2,x-^/2,z)e'rtdt;, (4) 

where xi = x + ^/2 andx2 =x — ^/2. The angle brackets (...) 
denotes statistical averaging which will be used when con- 
sidering extensions to partial coherence, see section [V] The 
Wigner transform method is frequently used in quantum the- 
ory and has also been successfully applied to the study of par- 
tially coherent light propagation in noninstantaneous nonlin- 
ear media. In the quasi-classical limit the formalism is further 
similar to the collision-less kinetic theory of plasma physics. 

Given the correlation or Wigner matrix it is easy to con- 
struct the Stokes vector S = S(x) which characterizes the po- 
larization state by its mapping onto the Poincare sphere. In- 
troducing the extended set of Pauli matrices Oj (<7o = I) we 
have 

Sj(x) = tr(r mn (x u x 2 )<Jj)\ xi=X2=x (5) 

for the correlation matrix, with tr(...) denoting the trace 
which is equivalent to 

Sj{x) = J tr{W mn <jj) dp (6) 

using the Wigner matrix. This agrees with the customary 
definition of the Stokes vector for circularly polarized light 
l26l l27tl . In component form we have 

SoMH|v/ + | 2 )(*) + <|v/-| 2 >(*) ( ? > 

Si(x) =2(v/+r-Cos(5))(jc) (8) 
S 2 (x)=2(v/+V/-sin(5))W (9) 

s 3 W = (|v/+l 2 )W-(|v/-l 2 )W, (10) 

with 8 = 0_ — + being the phase difference between the 
complex fields. 



III. SCALAR WIGNER FORMALISM 

To construct a scalar Wigner function we assume it to take 
the general form of a linear combination of the elements of 
the Wigner matrix. Specifically we exploit the analogy with 
spin- 1/2 systems in the Wigner formalism of quantum theory 
112 811 and define the scalar Wigner function as 

W(x,p,z,s) = Y<\ KM); • ^] m „W m „. (11) 

m.n 

Note that this choice of representation is invertible but not 
unique. Another choice for a scalar Wigner function describ- 
ing polarized light has previously been introduced by Luis in 
[12911 - The representation Eq.dTTb has the advantage of param- 
eterizing the unit Poincare sphere using the position vector 
s and maps the circularly polarized states onto the poles of 
the sphere. It should be emphasized that we do not treat the 
polarization as a distribution but rather take the vector § to 
determine the projection onto a particular state. 

Using the definition of the Wigner transform Eq.(|4|l it is 
easily found that the right- and left-handed circular intensities 
are given by 

/GO 
W(x,p,z,t)dp (12) 
-00 

and 

(\W-\ 2 ) = j~W{x,p,z-%)dp (13) 

while 

2^ W(x,p,z,x) dp = (|v/+| 2 ) + (\W-\ 2 ) +2(y/+y/--cos(<5)) 

(14) 

gives the intensity of the total field, including the beating term. 

The scalar Wigner function W(x,p,z,s) can be seen as de- 
scribing an instantaneous state of the system in an extended 
phase space. The evolution equation for this state is obtained 
by applying the Wigner transform to Eqs.(fl]i and (O, and then 
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using the definition Eq.dTTb- The particulars of similar deriva- 



tions have previously been presented elsewhere H30L 13111 . thus 
we will only present the result of the calculation which is 
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where we have introduced the difference between the propa- 
gation constants, A/3 — /3+ — /3_, and the intensities are pro- 
vided by the relations Eqs.dT2b and (fT3l . Note that the trigono- 
metric functions in Eq.dl3l should be interpreted as operators 
defined by their respective series expansions, with the arrows 
indicating the direction of application of the derivatives. 

Since the right- and left-handed circularly polarized inten- 
sities do not figure independently in Eq.dTSb but only in the 
combination of sum and difference intensities, it is convenient 
to write these in a shorthand form as 

/0 = <|V + | 2 } + (|V-| 2 } 

= 2/ W(x,p,z,0)dp = 2 I (l-s-Vg)W</p (16) 



and 



/3 = (|V+| 2 )-(|V-| 2 >=2 / z-VtWdp 
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The evolution equation Eq.(IT5b together with the relations for 
the sum and difference intensities thus constitute a closed sys- 
tem of partial integro-differential equations. Although this 
system appears quite complicated, it can directly furnish pro- 
jections onto a relevant evolution equation for arbitrary po- 
larization states s. Additionally, the formalism highlights the 
role played by unequal group velocities v g+ ^ v g _ and un- 
equal mode intensities I?, ^ 0, which may not be obvious from 
Eqs.([TJi and (|2). The system can in general be shown to pos- 
sess three invariants associated with the conservation of en- 
ergy, of momentum and of the Hamiltonian, viz. 
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These invariants are independent of the polarization direction 
and corresponds to the constant phase, x- and z-translation 
symmetries of the system Eqs.(Q]) and (|2]l by application of 
Noether's theorem. The evolution equation (TT~5T > is not known 
to have additional invariants, except in certain limits of the 
parameters, c.f. lf32ll . 

In the Vlasov limit the arguments of the trigonometric op- 
erators are small, i.e. (d/dx),(d/dp) <C 1. Keeping only the 
lowest order contribution in an expansion, it turns out that W 
is conserved along the orbits in the extended phase space. 

As a consequence Eq.(|T5b can be written as 
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(24) 

Since both d 2 f/dxdp = and Vg • % = 0, we may alterna- 
tively write Eq.d2TI) as a continuity equation in the extended 
phase space, viz. 

^^■(XW) + Tx (i p Wy Tp ^w)=0. (25) 

Which is an immediate consequence of the absence of photon 
creation and annihilation mechanisms. 



IV. MODULATIONAL INSTABILITY 

It is well known that stationary monochromatic solutions of 
the system Eqs.(Q]) and (f2]i can experience modulation insta- 
bility both in the anomalous (ay < 0) and normal (ay > 0) 



dispersion regimes. However, the stability analysis is com- 
plicated due to the general stationary continuous wave (CW) 
solution of Eqs.([T]i and (ffj being given in terms of Jacobian 
elliptic functions [25Q. Here, the Wigner transform method of- 
fers a convenient shorthand for obtaining a general dispersion 
relation, which illustrates many features of the polarization 
dependence, without the need to actually specify the back- 
ground solution. For simplicity we consider the equal group 



velocity case v g 



v g - throughout. 



To study the modulational instability associated with 
Eq.(|T5b, c.f. l33tl . we linearize the evolution equations of 
the Wigner matrix elements using the ansatz W m „ = W® m + 
£Wj,„ exp [i(kx — (Dz)} and consider the first order perturbation 
of the stationary CW background solution W® n as e — > 0. To 
cast the resulting dispersion relation in a form suitable for the 
chosen Wigner distribution function we make the simplifying 
assumption of low birefringence, i.e. Aj3 = 0. After some cal- 
culations it is then found that the dispersion relation can be 
written as 
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where the overlines indicates integration with respect to p, 

A/(p) = f(p + k/2) - f(p - k/2), Lf(p) = f(p + k/2) + 
f(p - k/2) and 77 = k (± + CCp) - CO. 

Although, the general CW solution of the system Eqs.(Q~|i 
and (O is given in terms of elliptic functions it is also possible 
to obtain solutions in terms of elementary functions when the 
light is polarized along either the fast or slow axis ll34ll . The 
solution of Eq.(fT3t then takes the following simple form 



W = I (lTs x )8{ P ), 



(27) 



which is seen to be independent of s y and s z . 

Another simplification is obtained for a purely incoherent 
interaction when the linear coupling term is neglected ( K = 0). 
The dispersion relation Eq.(l26b then reduces to 
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Considering now a general coherent CW background solution 



I 

of the form 
1 



W 







[(Iv+I 
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+ is y (yly + -y* +W _)+s z (\w + \ 2 -\y-\ 2 )]8(p). (29) 

It is seen that the modulational instability depends only on the 
total of the sum and difference intensities, regardless of the 
composition of the fields, since 



(l-8.V,)W = ijb5(p) J 



2-V 8 W = -/ 3 5(p), (30) 



where both Iq and ^3 are constant since the modes energies are 
now individually conserved. 

In this case the dispersion relation Eq.d28b can be evaluated 
explicitly with the result 



1 0*1 



(31) 

Eq.OTli clearly shown that the modulational instability can oc- 
cur also in the normal dispersion regime due to the nonlinear 
cross phase modulation, since for < < Iq we have that 
(1 + v) 2 / 2 > (1 + v) 2 / 2 -4v/ 3 2 > (1 - v) 2 / 2 . The instability 
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growth rate is given by 



again be evaluated in a closed form with the result 



\ak\ 2y 



a 



(l + v) 2 / 2 -4v/f-(l-v)/ 



(32) 

c.f. the corresponding expression for the partially coherent 
case in [33]. In the anomalous dispersion regime the instabil- 
ity growth rate analogously becomes 



|ce&| 2y 



a 



(1-v)/o+a/(1 + v) 2 / 2 -4v/ ; 



-k 2 . 



(33) 

In either case the maximum region of instability and the max- 
imum growth rate are obtained when the mode intensities are 
equal, i.e. /3 = 0. 



V. PARTIAL COHERENCE 



The above considerations can easily be extended to include 
the effects of partial coherence by allowing the fields to have 
a stochastically varying phase with statistical properties de- 
scribed by the mutual correlations between the fields 13511 . 
The partial coherence manifest itself as an additional width 
in p-space of the Wigner distribution function, beyond the in- 
trinsic spectral width of the coherent fields. In general, the 
coherence properties of each independent element of the her- 
mitian correlation matrix needs to be specified separately. 

When considering modulational instability, we assume, for 
analytical simplicity, that the statistical properties of the CW 
background solution for the two counter-rotating fields are 
identical, as would be the case if e.g. they were generated by 
the same source. We further assume that their coherence prop- 
erties can be described by a Lorentzian background distribu- 
tion, viz. 



(1-S-V S )W° 



P0 



2k p 2 
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2k p 2 
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■Pa 
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where p 1 is a parameter which can be identified with the 
correlation length. Note that in the coherent limit of po — >• 
0, the Lorentzian approaches a delta function and Eq.d30b is 
recovered as necessary. 
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For the case of K = 0, the dispersion relation Eq.(f28b may 



= 4r 2 + 1 (1 _ v)/o ± f\A l+v ) 2/ o- 4v/ 3-<>°' 

c.f. Eq.(|3"TTi and [33]. The partial coherence thus acts to sup- 
press the instability growth rate in a similar fashion to the 
Landau damping of electron plasma waves. Indeed, for a suf- 
ficiently large degree of incoherence the instability may be 
completely quenched. 

VI. CONCLUSIONS 



We have proposed a scalar Wigner distribution function 
for describing polarized light propagating in a nonlinear Kerr 
medium. By employing a circular polarization basis the polar- 
ization state is included in a manner similar to the spin state in 
the quantum kinetic theory of spin- 1/2 particles. The approach 
has been applied to optical propagation in a Kerr medium gov- 
erned by a coupled system of nonlinear Schrodinger equations 
and an extended phase space evolution equation has been de- 
rived along with invariant quantities. We have further con- 
sidered the problem of modulational instability and derived a 
dispersion relation valid for general background distributions, 
which in particular is shown to agree with the expected result 
in the limit of low birefringence and no linear coupling. Fi- 
nally the extension to partially coherent fields and its effect on 
the modulational instability has been considered. 

Future development of the above results could be in the di- 
rection of a 3D theory, as well as allowing for structurally 
different types of interactions (i.e., leaving the limitations of 
Kerr media). Such developments could prove of interest in 
e.g. magnetization studies through short laser pulses, where 
the statistical properties of the incoming light, and its inter- 
action with the matter target, can give rise to interesting new 
dynamics of spin structures. Thus, the coupling between dif- 
ferent degrees of freedom in light and matter would be able 
to capture a large set of fast dynamics. Moreover, such gen- 
eralizations could also include a more generic description of 
the invariants of the field. Such a development would allow 
for a less than ad hoc formulation of the spin and angular mo- 
mentum properties of light in terms of a scalar distribution 
function. It is well known that the angular moment of light 
can readily affect the dynamics of matter structures; it would 
be of interest to see how such a formulation could highlight 
the interaction between partial coherence, spin/angular mo- 
mentum, and matter dynamics induced by light. 
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